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NOMENCLATURE 


pulse  frequency,  pulses /sec 

transfer  function  of  component  in  control  loop, 
dimensionless 

2 

moment  of  inertia  (e.  g.  ,  about  pitch  axis),  slug-ft 

2 

moment  of  inertia  about  roll  axis,  slug-ft 

2 

moment  of  inertia  about  yaw  axis,  slug-ft 

thrust  impulse  single  corrective  torque  application,  lb- sec 

total  thrust  impulse  over  operating  time,  T,j,,  lb-sec 

stabilizing  or  destabilizing  torque  parameter,  dimensionless 

position  gain,  dimensionless 

rate  gain,  sec 

gain  of  i*^  component  in  control  loop,  dimensionless 
moment  arm,  ft 

modulation  factor  (m  =  1  -*  constant  thrust),  dimensionless 
threshold  modulation  factor,  dimensionless 
steady  external  torque,  ft -lb 

stabilizing  or  destabilizing  external  torque  coefficient, 
ft -lb/ radian 

number  of  pulses  required  to  reverse  angular  rate, 
dimensionless 

total  number  of  thrust  impulses  in  time,  T,^,  dimensionless 
rate  measurement  threshold,  radians/sec 
Laplace  operator,  sec~^ 
time,  sec 
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thrust  pulse  width,  sec 

minimum  thrust  pulse  width,  sec 

coast  time  w’ith  destabilizing  external  torque,  sec 

coast  time  without  external  torque,  sec 

coast  time  with  destabilizing  external  torque,  sec 

coast  time  with  steady  external  torque,  sec 

thrust  pulse  magnitude,  lb 

system  operating  time,  sec 

thrust  generator  threshold,  dimensionless 

total  propellant  used  in  time  T^,  lb 

position  measurement  threshold,  radians 

error  signal,  radians 

vehicle  Euler  angle,  radians 

position  measurement  signal,  radians 

Laplace  transformed  position,  radians 

initial  condition  in  0  (beginning  of  coast),  radians 

position  command  signal,  radians 

vehicle  angular  rate,  radians/ sec 

rate  measurement  signal,  radians/sec 

initial  condition  in  0  (beginning  of  coast),  radians/sec 

angular  rate  reversal  magnitude,  radians/ sec 

_2 

stabilizing  or  destabilizing  torque  parameter,  sec 
steady  external  torque  parameter,  radians/ sec^ 
critical  value  of  steady  torque  parameter,  radians/sec 


IV 


relative  propellant  consumption  -  destabilizing  torque, 
dimensionless 

relative  propellant  consumption  -  destabilizing  torque, 
dimensionless 


relative  propellant  consumption  -  stabilizing  torque, 
dimensionless 


relative  propellant  consumption  -  steady  torque, 
dimensionless 

orbit  frequency,  radians /sec 


V 


I.  DESCRIPTION 


The  single -axis  attitude  control  system  described  by  the  diagram  in 
Figure  1  is  now  considered.  The  position  and  rate  measurements  are  con¬ 
sidered  to  have  independent  and  symmetric  threshold  values  6  and  r,  re¬ 
spectively.  A  threshold  error  signal  of  magnitude  v  is  necessary  to  obtain 
an  output  from  the  torque  generating  device. 

The  authors  consider  in  particular  that  the  thrust  which  generates 
control  torque  is  pulse  modulated  and  the  modulation  factor  is 


m  =  fAt 


It  is  possible  to  modulate  the  thrust  with  pulse  frequency  f  or  pulse 
width  At,  or  both.  Figure  2  depicts  the  meaning  of  these  quantities  along 
with  T^  the  thrust  pulse  magnitude.  The  time  average  thrust  will  obviously 
saturate  for  m  =  1  but  this  condition  will  not  be  a  point  of  concern  in  this 
paper. 
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Fig.  1.  Attitude  control  sy.<»tem  with  nonlinearities. 


II.  SWITCH  LINES 


The  following  relations  can  be  written  from  Figure  1  by  inspection 


10,  jej  <  5  (la) 

0-6,  0  >  6  (lb) 

0  +  8,  0  <  -  6  (Ic) 

10,  I©!  <  f  (2a) 

0  -  f,  0  >  f  (2b) 

0+f,0<-r  (2c) 

10,  |« I  <  V  (3a) 

Kd  (<  -  v)  +  m^,  *  >  V  (3b) 

Kjj  (c  +  v)  -  m^,  e  <  -  V  (3c) 


-t 


0 


+ 


0 


(4) 


The  equations  of  the  lines,  referred  to  as  switch  lines,  in  the  phase 
plane  that  separates  the  coasting  region  from  regions  where  corrective  torque 
is  applied  may  be  determined  by  setting  m  ■±  m^  =  0  in  Eqs.  (3b)  and  (3c). 
The  equations  of  the  switch  lines  can  best  be  found  by  constructing  a  table  of 
all  the  combinations  of  signs  of  0,  0,  and  c  (see  Table  I).  A  set  of  equations 
for  the  0  and  0  signals  at  or  above  threshold  are  obtained  from  Eqs.  (lb), 
(Ic),  (2b),  (Zc),  (3b),  (3c),  and  (4)  (see  column  B,  Table  I).  Another  set  is 
obtained  by  letting  0  drop  below  threshold  or  (0  ±  r)  =  0  in  the  equations  of 
column  B  (see  column  C).  The  remaining  possible  equations  are  obtained  by 
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Table  I.  Equations  defining  switch  lines  in  the  phase  plane. 


0 


4 


letting  6  drop  below  the  threshold  or(0±6)=Oin  the  equations  of  column 
B  (see  column  D).  Some  of  the  combinations  of  signs  are  not  possible  and 
are  so  indicated.  This  happens  because  of  the  inequalities  in  the  mother  set 
of  equations.  The  segments  of  the  lines  defined  by  the  equations  in  columns 
B,  C,  and  D  w'hich  enclose  the  coast  region  are  plotted  in  the  phase  plane 
diagram  of  Figure  3. 
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III.  POWERED  TRAJECTORIES 

(7 


It  is  seen  that  the  switch  lines  defined  by  Figure  3  and  Table  I  separate 
regions  in  the  phase  plane  where  the  trajectories  differ  in  character.  Outside 
the  coast  region  are  powered  regions  where  the  trajectories  are  those  where 
corrective  torqvie  is  applied.  The  following  trajectory  descriptions  apply  to 
the  system  under  consideration  without  external  torques. 

Regions  1,  7,  8,  and  14  contain  rectilinear  portions  of  the  trajectories 
where  the  control  system  operates  without  benefit  of  position  feedback.  Re¬ 
gions  2,  5,  6,  9.  12,  and  13  contain  the  second  order  system  spiral  portions 
of  the  trajectories.  Regions  3,  4,  10,  and  11  contain  the  elliptical  portions 
of  the  trajectories  where  the  system  operates  without  benefit  of  rate  feedback. 

The  regions  of  importance  so  far  as  limit  cycles  are  concerned  are  the 
coast  region  and  the  elliptical  trajectories  regions,  for  it  can  be  deduced 
that  the  steady  limit  cycle  of  the  system  being  considered  is  determined  by 
f  and  (6  +  v).  The  vehicle  will  coast  at  angular  velocity  ±  f  until 
0  =  ±  (6  +  v),  at  which  time  corrective  torque  is  applied  without  benefit  of 
rate  feedback,  giving  rise  to  an  elliptical  trajectory  symmetrical  about  the 
0  axis.  The  corrective  torque  is  removed  when  the  point  on  the  elliptical 
t  rajectory  coincides  with  a  point  on  a  switch  line.  It  should  be  pointed  out 
that  the  powered  portion  of  the  limit  cycle  trajectories  is  elliptical  only  if 
the  threshold  modulation  factor  m  =0.  If  m  ^0,  the  trajectories  tend 
to  be  parabolic  and  are  perfect  parabolas  if  m  =  m^  =  constant.  When  pulse 
modulation  is  used,  the  minimum  angular  momentum  impulse  applied  to  the 
vehicle  by  a  single  thrust  pulse  should  be  much  less  than  If  or 


T  iAt.  =IAe<<If 
m  min 


(5) 
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in  order  that  the  foregoing  comments  apply.  If  this  is  not  the  case,  the 
mean  drift  velocity  will  be  determined  by 

T  iAt  . 

A  n  m  min 

e  =  2  -  — -  (6) 

where  n  is  the  integer  number  of  pulses  required  to  reverse  angular  rate. 
Since  n  should  be  an  integer,  it  will  be  some  value  higher  than  the  reversing 
impulse  or 


"iT-Tsr—  <’> 

m  min 

Ideally  n  should  be  the  next  integer  value  above  the  value  of  the  right  side 
of  the  above  expression.  An  interesting  thing  to  note  is  that  the  number  of 
pulses  should  lie  between  that  required  to  insure  a  change  of  sign  of  0  and 
that  required  to  change  6  from  f  to  -f.  Thus 


T - TET—  <  ^  ^  T - iAt  . 

m  min  m  min 

and,  therefore,  the  mean  drift  velocity  will  be  less  than  f.  For  n  large, 
the  torque  control  becomes  almost  continuous  (small  granularity),  n  ap¬ 
proaches  the  right  side  of  Eq.  (8)  and  the  drift  velocity  is  approximately  f. 
For  n  =  1  and  Ir/T  f  At  nearly  equal  to  but  less  than  unity,  the  mean 
drift  velocity  approaches  r/2.  It  would  appear  that  by  proper  design,  an 
attitude  control  system  using  pulse  modulation  might  have  a  built-in  safety 
margin  of  gas  or  propellant  supply  (in  the  absence  of  external  torque). 
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Figure  4  illustrates  a  pulse  frequency  modulation  system  response 
trajectory  in  the  phase  plane  where  the  mean  drift  velocity  is  slightly  greater 
than  f/2.  The  pulse  width  in  this  case  is  constant.  Figure  5  illustrates  a 
pulse  m-odulation  system  response  trajectory  when  the  pulse  width  and  pulse 
frequency  are  both  modulated.  Position  and  rate  thresholds  are  not  shown 
in  this  figure. 
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Fig.  4.  Pulse  frequency  modulation  system  response. 
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IV.  COAST  TRAJECTORIES 


In  the  absence  of  external  torques,  the  vehicle  will  coast  with  constant 
angular  velocity  across  the  coast  region  from  switch  line  to  switch  line  and 
from  left  to  right  in  the  upper  half  of  the  phase  plane  and  from  right  to  left 
in  the  lower  naif  of  the  phase  plane. 

The  characteristics  of  the  coast  trajectories  when  the  vehicle  is  subject 
to  external  torques  can  be  determined  by  solving  the  vehicle's  equation  of 
motion.  The  equation  of  motion  of  the  coasting  vehicle  with  a  stabilizing  or 
destabilizing  torque  and  a  steady  torque  acting  on  it  is 


I  6 


M 

s 


(9) 


where  8M/30  =  stabilizing  or  destabilizing  torque  coefficient  (ft-lb/radian) 
and  =  steady  torque  (ft/lb).  The  coefficient  9M/30  is  defined  as  the 
moment  (stabilizing  or  destabilizing)  per  unit  angle  0. 

Considering  aerodynamic  and  gravitational  moments  acting  in  the  pitch 
plane  on  a  circular  orbital  vehicle  witn  pitch  attitude  deviation  0  from  the 
local  horizontal  in  the  direction  of  motion,  the  torque  coefficient,  3M/90, 
can  be  expressed  as 


3M 

w  " 


^0 


(I  -  I  ) 

O  X  z 


where  Mg  is  the  aerodynamic  stability  derivative  with  respect  to  0,  I^  is 

the  yaw  moment  of  inertia,  I  is  the  roll  moment  of  inertia,  and  w  is  the 

X  o 

vehicle  orbit  frequency.  The  second  term  on  the  right  of  this  relation  is  the 
torque  per  unit  angle  0  resulting  from  the  gravity  gradient  as  presented  and 
discussed  by  Frye  and  Stearns  (Ref.  1).  As  is  pointed  out  by  Frye  and 
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Stearns,  Davis  (Ref.  7)  has  shown  that  this  term  arises  from  a  combination 
of  the  effects  of  torque  due  to  gravity  gradient  and  centrifugal  force  of 
orbital  motion.  More  general  relations  for  the  gravity  torques,  from  which 
one  can  determine  the  gravitational  torque  derivative,  are  available  in  the 
literature  (Refs.  3,  4,  8,  13,  14).  The  aerodynamic  stability  derivative 
may  be  deduced  from  the  normal  force  and  axial  force  characteristics  of 
the  vehicle  resulting  from  what  may  be  considered  Newtonian  or  free  mo¬ 
lecular  flow  as  is  discussed  by  Frye  and  Stearns  (Ref.  1),  and  DeBra  and 
Stearns  (Refs.  5,  6). 

A  steady  torque  with  respect  to  body  axes  could  arise  from  solar 
radiation  pressure  or  winds.  In  this  case,  the  torque  would  be  "steady" 
only  during  the  time  of  year  when  the  orbit  plane  of  an  earth  satellite  would 
be  normal  to  Earth-sun  radius  vector.  The  solar  pressure  torque  would  be 
constant  for  a  vehicle  in  a  circular  heliocentric  orbit  where  one  axis  is  kept 
aligned  with  the  velocity  vector.  It  can  also  arise  from  gas  ventage,  leakage, 
and  from  steady  aerodynamic  moments  in  a  vehicle  whose  axis  of  symmetry 
is  not  nominally  aligned  with  the  orbital  velocity  vector.  (Refs.  1,  4,  10, 

11,  12) 

Equation  (9)  may  be  written  as 


e 


0  +  p 


(10) 


where 


2  A  8M/8e 

I 


(11) 
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If  Eq.  (10)  is  divided  by  6 


0 

6 


de 


+  ti 


0 


(12) 


Equation  (12)  integrates  directly  to 

2  2  2 

0  -  X  0  -  2|i0  =  constant 


(13) 


2 

If  the  sign  of  X  is  positive  and  p  =  0,  the  external  torque  is  destabilizing 
and  the  result  is  the  equation  of  a  family  of  hyperbolas  with  center  at  the 
origin  and  asymptotes  with  slope  ±  X.  These  hyperbolic  trajectories  exist 
only  within  the  coast  region  of  the  phase  plane.  Figures  6  and  7  show  the 
two  types  of  limit  cycles  which  are  possible  assuming  that  higher  order 
effects  will  preclude  steady  limit  cycles  inside  the  ones  shown.  In  Figure  6 
the  coast  trajectories  resemble  the  constant  0  condition  which  is  obtained 
with  no  torque.  In  Figure  7  the  coast  trajectories  are  such  that  0  always 
h'as  a  mean  value  other  than  zero  and  never  reaches  zero  (unless  the  trajec¬ 
tory  follows  the  asymptotes). 

2 

If  the  sign  of  X  is  negative  and  p  =  0,  the  external  torque  is  stabi¬ 
lizing  and  we  have  the  equation  of  a  family  of  ellipses  which  describe  the 
trajectories  in  the  coast  region.  Figure  8  illustrates  the  lowest  frequency 
steady  limit  cycle  possible  with  stabilizing  torque.  Elliptical  trajectories 
are  possible  inside  the  one  shown  but  it  is  assumed  that  spurious  torques 
and  lags  in  a  practical  system  will  eventually  open  these  out  to  the  one  shown. 

2 

If  X  =0  and  p  is  positive,  the  equation  of  a  family  of  parabolas  is 
obtained.  These  parabolas  open  to  the  right  in  the  phase  plane.  Figure  9 
illustrates  the  steady  limit  cycle  obtained  assuming  that  higher  order  effects 
will  open  out  the  trajectories  that  are  possible  inside  the  one  shown.  For  p 
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smaller  than  that  of  Figure  9,  the  coast  trajectory  flattens  until  the  apex  of 
the  parabola  coincides  with  the  9  =  -  (6  +  v)  switch  line  and  further  re¬ 
duction  of  [j.  causes  corrective  torque  application  on  the  left  side  of  the 
phase  plane. 

Had  the  steady  torque  term  p  not  b^en  neglected,  in  the  destabilizing 
and  stabilizing  external  torque  types  of  ♦^’•aiectories  mentioned  in  the  fore¬ 
going,  its  effect  would  have  been  to  shift  the  hyperbola's  and  ellipse's  centers 
along  the  0  =  0  axis.  Figures  lu  and  il  show  this  effect  and  the  resulting 
limit  cycle  trajectory  shapes. 
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ig.  10.  Destabilizing  torque  plus  steady  p-g  ^  Stabilizing  torque  plus  steady 

torque  trajectory.  torque  trajectory. 


V.  COAST  TIME 


The  ratio  of  the  time  of  corrective  torque  application  to  the  drift  or 
coasting  time  in  a  practically  designed  system  will  likely  be  of  the  order  of 
10  percent  or  less.  Thus  the  period  of  the  limit  cycle  in  most  cases  will  be 
dependent  primarily  upon  the  drift  time.  This  time  can  be  determined  for  a 
system  with  external  torques  and  can  be  compared  with  the  drift  time  of  a 
system  without  external  torques  to  obtain  an  estimate  of  relative  propellant 
consumption  since  the  amount  of  propellant  expended  over  a  time  interval  of 
operation  of  the  system  will  be  proportional  to  the  frequency  of  the  limit 
cycle. 

Equation  (10)  may  be  solved  subject  to  the  initial  conditions  that  exist 
at  the  beginning  of  the  drifting  portion  of  the  limit  cycle. 

0  (o)  =  e 

o 

and  (14) 


e  (o)  =  0 
'  '  o 


Definition  is  given  of 


L  |0(t)|  =  §(8) 


Taking  the  Laplace  transform  of  Eq.  (10) 


0 


§  .  . 

0- 

o 


8  + 


— 1 - ^  ^ - T 

{3^  -  \^)  0^8  (8^  -  \^) 


(15) 
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Finding  the  inverse  Laplace  transform  of  Eq.  (15)  obtains  for  X,  positive 


o 


(t)  =  cosh  \  t  +  sinh  X  t  +  —  [cosh  X  t  -  l]  ,  X  >0  (16) 


~T 
\  Q 


and  for  X  negative 


w 

cos  |X|  t  +  sin  |X|t  +  j-|j| —  [  1 

II  o  |X  0 

•  r\ 


-  cos  |X|t]  ,  X  <  0  (17) 


The  coast  time  interval  from  the  initial  conditions  to  6  =  0  is  obtained  by 
setting  0  (t)  =  0. 

2 

A.  Destabilizing  External  Torque,  X  >0  with  p  =  0 

The  coast  time  between  corrective  torque  application  from  Eq.  (16) 
for  the  limit  cycle  type  shown  in  Figure  6,  where  p  =  0,  is 


*  2  ^  ,  -1 
'Dj  =  \ 


(18) 


where 


K 


(19) 


Using  the  identity 


tanh 
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Eq.  (18)  can  be  expressed  as 


(20) 


which  corresponds  to  the  condition  K  >  1. 

The  coast  time  between  corrective  torque  application  for  the  limit 
cycle  shown  in  Figure  7  is  determined  by  differentiating  Eq.  (16)  and  setting 


obtaining 


t 


D 


2 


(21) 


Which  corresponds  to  the  condition  0  <  K  <  1.  . 

Comparing  Eqs.  (20)  and  (21)  it  can  be  seen  that  in  general 

0 


D 


=  iin 


IJS. 


+  1 


iK  -  1 


K  > 


(22) 
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B.  Stabilizing  External  Torque,  X  <0  With  p  =  0 

2 

The  coast  time  between  the  switch  lines  for  X.  <0  and  p  =  0  is 
obtained  from  Eq.  (17) 


tan 


-1 


(23) 
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2 

C.  Steady  External  Torque,  |jl  ^  0  with  X.  =  0 

The  coast  time  for  the  type  of  trajectory  shown  in  Figure  9  can  be 
found  by  differentiating  Eq.  (17)  and  setting  X  =  0  and  0  =  0.  This  results 
in  the  drift  time  between  corrective  torque  applications 


(24) 


When  p  is  small  enough  that  corrective  torque  is  applied  on  both  sides 
of  the  phase  plane  diagram,  Eq.  (17)  is  again  differentiated,  setting  X  =  0 
and  setting 


0  = 


(25) 


to  find  the  drift  time.  Equation  (25)  is  obtained  from  Eq.  (13)  with  0  =  -0^ 
where  the  constant  is  evaluated  using  the  initial  conditions  0=0^  and 

0  =  -0  .  The  result  is  obtained 

.  o 


I 


t  = 


(26) 


for  the  coast  time  between  corrective  torque  applications.  The  value  of  p 
for  which  the  apex  of  the  parabolas  coincides  with  G  =  -  0^  is  obtained 
from  Eq.  (25)  with  0  =  0.  It  is  what  may  be  termed  the  critical  value  of 
the  steady  torque  parameter 


(27) 
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VI.  RELATIVE  PROPELLANT  CONSUMPTION 


If  it  is  assumed  that  corrective  torque  application  time  is  negligible 
compared  to  coast  time,  the  propellant  consumption  will  be  proportional  to 
corrective  torque  application  frequency  or  inversely  proportional  to  coast 
time  for  a  given  operating  time.  The  relative  propellant  consumption  0  is 
defined  as  (propellant  consumption  with  external  torque)  +  (propellant  con¬ 
sumption  without  external  torque). 

A.  Destabilizing  External  Torque 

For  the  case  of  destabilizing  external  torque 


A 


From  Eq.  (21),' using  L'Hospital's  rule 


lim 

X-0 


^t 


(28) 


(29) 


which  could  have  been  deduced  directly. 

Using  Eqs.  (22)  and  (29)  in  Eq.  (28) 

»D = ij t^n i »  <'»> 

It  is  of  interest  to  see  how  propellant  consumption  varies  with  the 

2 

destabilizing  torque  coefficient.  Accordingly  ^  is  plotted  vs.  1/K  in 
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Figure  12.  This  figure  shows  a  theoretical  relative  propellant  consunaption 
of  zero  for  K  =  1.  Practically,  less  impulse  than  for  the  case  of  a  system 
without  torque  might  be  expected  if 

0  <  -^  <  1.  44 


or 


(31) 


0.  83  <  K  <  oo 

However,  this  presumes  that  the  point  (6^.  6^)  which  corrective 
torque  is  applied,  is  the  same  for  both  conditions  (i.e.  ,  with  torque  and 
with  no  torque).  This  behavior  is  seen  in  the  dashed  curve  of  Figure  13. 

The  case  is  examined  where  there  is  a  given  value  of  torque  and 
hence  It  should  be  possible  to  vary  K  by  virtue  of  the  possibility  of 
adjusting  The  question  may  reasonably  be  asked:  What  effect  does 

the  choice  of  switching  point  {6q»  9^)  have  on  propellant  consumption?  To 
answer  this  it  may  be  noticed  from  Figure  12  that  K  =  0.  83  gives  the  same 
propellant  consumption  as  the  case  of  zero  torque  (it  may  be  considered 
unity  consumption).  The  ratio  is  formed 


,  .  -0/0  =  0.  83\ 

'A  o  o 

=  ■ 


f  n 


K  +  1 


K  -  1 


K  =  0.  83 


0  =  K\ 

o  o 


k  IK 


K  +  1 


(32) 
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Q 


or 


I 


2.  4 


in 


1 

WTI 

K  -  1 


(33) 


This  relation  is  plotted  as  curve  1  of  Figure  13  and  is 


propellant  consumption  when  switch  "id  t 
propellant  consumption  when  switched  at  =  0-  83X, 


propellant  consumption  when  switched  at 
propellant  consumption  when  switched  at  with  X  =  0 


From  the  solid  curve  of  Figure  13,  one  sees  that  for  a  given  value  of 

torque 


0.  83 


1 

"  T7I 


<  K  <  1.  2 


(35) 


in  order  that  consumption  be  less  than  that  obtained  with  zero  torque.  It  is 
apparent  from  Eq.  (28)  that 


^  (K)  = 


(36) 


Something  might  be  gained  from  this  observation;  i.  e.  ,  if  the  torque  coeffi¬ 
cient  is  predictable  within  a  ±  percent  expected  error,  it  would  be  better  to 
set  K  >  1  because  the  penalty  for  error  is  not  so  great  in  this  direction  as 
for  K  <  1. 


"-I 
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In  regard  to  the  means  for  adjusting  one  observes  from  Fig 

ure  3  that  6  =  r  and  6  =  6  +  v  and,  therefore 

o  o 


-  TTT 

o 


(37) 


As  a  result,  the  switching  point  is  easily  adjusted  if  these  quantities  are  at 
one's  disposal  in  a  practical  system. 

B.  Stabilizing  External  Torque 

2 

The  relative  propellant  consumption  for  the  case  where  X  is  negative 
and  |i  =  0  is  obtained  using  Eq.  (23) 


t  „  ^ 
s  K  tan 


1 

rr 


Wi 


(38) 


Equation  (38)  is  plotted  in  Figure  12  which  shows  that  the  relative  propellant 
consumption  increases  monotonically  and  the  consumption  with  a  negative 
stabilizing  external  torque  coefficient  will  always  be  greater  than  the  con¬ 
sumption  without  torque. 

C.  Steady  Torque 

2 

The  relative  propellant  consumption  for  the  case  where  X  =0  and 
is  obtained  using  Eq.  (24) 


(39). 
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To  obtain  the  relative  propellant  consumption  for  \  =0  and  p  positive 

and  less  than  it  is  assumed  that  for  the  corrective  torque  on  the  left  side 
of  the  phase  plane,  one  consumes  only 


as  much  fuel  in  reversing  0  as  on  the  right  side.  The  ratio 


I 


will  then  be  the  relative  fuel  consumption  of  this  cycle  due  to  the  unsymmetric 
limit  cycle.  Then  the  product  of  this  quantity  and  the  ratio  of  the  coast  time 


should  give  a  reasonable  approximation  to  the  relative  fuel  consumption.  Thus 
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Equations  (39)  and  (40)  are  plotted  in  Figure  14,  where  it  is  shown  that  the 
fuel  consumption  of  the  system  can  be  as  small  as  1/4  the  fuel  consumption 
without  external  torques. 
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STEADY  TORQUE  RELATIVE  PROPELLANT  CONSUMPTI 


Fig,  14.  Steady  torque  relative  propellant  consumption. 
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VII.  TOTAL  PROPELLANT  CONSUMPTION 
WITHOUT  EXTERNAL  TORQUE 


The  following  assumptions  are  made: 

1.  Steady  limit  cycle  only  contributes  to  propellant  consumption. 

2.  No  external  torques. 

3.  Corrective  torque  application  time  negligible  compared  to  coast 
time. 

4.  Reaction  jets  only  are  used. 

One  then  has  the  set  of  equations  for  a  single  axis 


t  = 
o 


e 

o 


^T  =  Nto 
=  Ntji: 


lAe 


=  T  At 
m 


Ae  = 
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Simple  substitutions  yield 


W 


T 


T 

o 


(41) 
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VIII.  HIGHER  ORDER  EFFECTS 


In  practical  systems  the  dynamics  of  the  hardware  used  in  the  control 
system  lead  to  higher  than  second  order  system  response.  Figure  15  rep¬ 
resents  a  typical  single  axis  system  with  the  transfer  function 


e 

T 

c 


G2  (8)  G3  (S) 


t)  4  ®  ^  °1 


1  + 


(s)  G^  (8)  G3  (s) 


Figure  16  presents  a  comparison  between  a  second  order  system, 
i.  e.  ,  Gj  (s)  G2  (8)  G3  (s)  equals  one,  and  a  higher  order  system.  The 
phase  plane  is  the  result  of  an  analog  computer  study  of  a  system  without  a 
threshold  in  the  position  channel  to  the  valve.  The  dynamics  of  the  control 
system  hardware  cause  the  time  of  switching  to  lag  that  of  an  ideal  second- 
order  system  thereby  increasing  limit  cycle  rates  and,  subsequently,  fuel 
consumption. 
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IX.  DESIGN  CONSIDERATIONS 


A.  Minimum  Pulse  Width  •• 

The  number  of  pulses  to  reverse  the  direction  of  motion  when  corrective 
torque  is  applied  needs  to  be  one  or  more.  Some  advantage  accrues  from  a 
single  thrust  pulse  as  was  discussed  in  the  section  on  Powered  Trajectories. 

If  N  =  1  in  Eq.  (8)  then 


Ir  <  T  fAt  .  <  2If 

m  min 


(42) 


A  pulse  width  can  be  found  which  will  cause  0  to  change  its  sign  and 
also  which  will  provide  some  margin  of  safety  suc^i  that  two  corrective  pulses 
will  not  occur.  Choosing  a  value  of  torque  impulse  midway  in  the  range 
[  Eq.  (42)]  one  could  use 


At 


min 


1.  5  If 

^ - j 


m 


(43) 


The  minimum  or  threshold  pulse  width  should  be  at  least  this  small. 
B.  Maximum  and  Minimum  Pulse  Frequency 

If  pulse  frequency  modulation  alone  is  used  it  is  noted  that 


f 

max 


1 


At  . 


min 


(44) 


when  m  =  1  (or  continuous  thrust).  The  minimum  pulse  frequency  should 

be  such  that  1/f  .  or  the  time  between  pulses  at  threshold  conditions  is 
min  ^ 

long  compared  with  the  instrument  (and  possible  shaping  network)  response 
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times.  This  is  desirable  to  prevent  application  of  more  corrective  thrust 
impulses  than  is  necessary  because  of  slow  information  processing  by 
control  system  components.  If  both  of  the  desired  conditions  (i.  e.  , 
and  cannot  be  met  because  of  frequency  range  limitations  of  pulse 

frequency  generating  devices,  it  may  be  desirable  to  use  a  combined  pulse 
frequency  and  pulse  width  modulation  method  in  order  to  extend  the  modu¬ 
lation  factor  range  of  the  pulse  producing  equipment. 

C.  Adjustments  for  External  Torque 

The  minimum  pulse  width  should  likely  be  reduced  when  the  vehicle 
has  external  torque  acting  on  it.  This  can  be  deduced  from  curve  1  of 
Figure  13;  i.  e.  ,  it  is  not  desirable  to  allow  the  single  torque  impulse  to 
govern  A0  (unless  it  corresponds  to  K  =  1),  It  is  probably  better  from  a 
design  point  of  view  to  require  something  of  the  order  of  10  pulses  to  achieve 
A6  =  2f  in  order  that  the  adjustment  factor  f /6  +  v  predominate  in  the 
determination  of  relative  propellant  consumption. 

Then  again,  if  f  is  not  accurately  known  or  is  time  varying  in  some 
random  way,  it  may  be  better  to  choose  to  correspond  to  the  optimum 

value  of  K  =  1,  providing  that  minimum  pulse  width  can  be  accurately  set 
and  held  time  invariant  better  than  f. 


35 


X.  CONCLUSIONS 


1.  In  the  threshold  region  or  in  the  vicinity  of  the  origin  in  the  phase 
plane,  trajectories  are  or  can  be  piecewise  linear,  hyperbolic, 
parabolic,  and  elliptic. 

2.  Minimum  pulse  size  can  affect  the  mean  drift  rate  advantageously  or 
disadvantageously. 

3.  There  may  be  an  advantage  in  modulating  both  pulse  width  and  pulse 
frequency  in  order  to  reduce  the  single  pulse  impulse  below  If /i  so 
that  it  is  not  the  controlling  factor  in  the  resulting  mean  drift  rate. 

4.  With  destabilizing  torque,  the  vehicle  may  drift  (hyperbolic  trajectories) 
from  positive  to  negative  attitudes  or  the  limit  cycle  may  be  such  as  to 
yield  only  positive  or  negative  values  of  attitude. 

5.  Propellant  consumption  with  destabilizing  torque  is  less  than  propellant 
consumption  without  external  torque  for  values  of 

6.  Propellant  consumption  with  destabilizing  torque  is  less  than  propellant 
consumption  without  external  torque  for  values  of 
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if  6  /0  for  the  case  of  no  external  torque  is 
o  o 


0.  83 


Propellant  consumption  with  stabilizing  torque  will  be  greater  than 
propellant  consumption  without  external  torque  for  a  steady  powered 
limit  cycle  with  amplitude  determined  by  the  threshold  magnitudes. 

Propellant  consumption  with  steady  torque  is  less  than  propellant 
consumption  without  external  torque  for 


o 


Higher  order  effects  or  additional  lags  in  a  single  axis  attitude  control 
system  cause  propellant  consumption  to  be  greater  than  a  simple 
second  order  system. 

Total  propellant  consumption  due  to  limit  cycles  in  the  system  without 
torque  is  proportional  to  the  square  of  the  rate  measurement  threshold 
and  inversely  proportional  to  the  position  measurement  threshold  plus 
valve  threshold. 

Minimum  pulse  width  should  be  made  compatible  with  the  rate  meas¬ 
urement  threshold. 

Minimum  pulse  frequency  should  be  low  compared  to  sensor  cutoff 
frequencies. 

Something  of  the  order  of  10  pulses  should  be  designed  for  in  the  time 
duration  of  corrective  torque  application.  This  should  enable  the  dead 
zone  settings  of  f ,  6,  and  v  to  independently  determine  the  character 
of  the  limit  cycles. 
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